Saddle-point van Hove singularity and the phase diagram of high-T c cuprates 
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We examine the generic phase behavior of high-T c cuprate superconductors in terms a universal 
van Hove singularity in the strongly overdoped region. Using a rigid ARPES-derived dispersion we 
solve the BCS gap equation and show that the pairing interaction or pairing energy cutoff must be a 
rapidly declining function of doping. This result is prejudicial to a phonon-based pairing interaction 
and more consistent with a magnetic or magnetically enhanced interaction. 

PACS numbers: 74.25.Dw, 74.25. Jb, 74.62.Dh, 74.72.-h 
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The high-r c cuprates are remarkable in that, despite 
their varied atomic structures and consequent variation 
in bare electronic band structure [lj, they exhibit univer- 
sal phase behaviour. For example, with the exception 
of La2- a: Sr a; Cu04 which seems more prone to stripe in- 
stabilities, the thermoelectric power (TEP) is a universal 
function of hole concentration, Moreover, the evolu- 
tion with doping of both the thermodynamic properties 
and the spin susceptibility is likewise universal Most 
importantly, the overall temperature-doping phase dia- 
gram seems to be universal [J]. The common features in- 
clude: the onset of superconductivity at p « 0.05, the lo- 
cation of so-called u l/8 th anomalies" &tp « 0.12, optimal 
doping at p w 0.16, critical doping where the pseudogap 
closes at p « 0.19, the pseudogap line T*(p), and possi- 
bly the superconductor/metal transition at p « 0.27. To 
these we now wish to add another apparently common 
feature, namely the presence of a van Hove singularity 
(vHs) in the heavily overdoped region 

The E(k) dispersion for the hole-doped cuprates ex- 
hibits a saddle-point singularity, the so called van Hove 
singularity, sited at (tt,0) on the Brillouin zone bound- 
ary. Increased hole doping moves the Fermi energy, Ep, 
down towards the vHs which it eventually crosses [5j and 
where the density of states (DOS) diverges. Within a 
BCS picture the transition temperature is given by[(| 



ksTc = 1.14furfp exp(- 
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so that if the vHs crossing occurs within the supercon- 
ducting domain the exponential dependence upon the 
DOS should result in a local peak in T c , precisely at the 
vHs. Of course the cuprates do exhibit a SC phase curve 
that passes through a peak at optimal doping[4[ and there 
have been many attempts to explain this phase behav- 
ior in terms of a vHs crossing there [7j, |S, l9j|. Further, 
the rise in isotope effect exponent with underdoping was 
also explored as a consequence of a proximate vHsflol]. 
Finally, underdoped cuprates exhibit a pseudogap in the 
normal-state DOS that, amongst other things, causes a 
strong suppression of the spin susceptibility, % s , at low 



temperature [1, 11 1. Some groups sought to explain this 
suppression in terms of a nearby vHs[12J where the Fermi 
window, for T > 0, extends to the far side of the vHs 
thereby reducing \s- 

These ideas all failed for various reasons. The value of 
Xs near a vHs never falls more than 10% at low tempera- 
ture and so the model could never account for the nearly 
full suppression of \s as T -c 00, [Tl|. Further, ARPES 
revealed that the vHs at (ir, 0) lay at least 60meV be- 
low Ep at optimal doping and thus could not account 
directly for the location of the maximum in T c . 

However, we now need to revisit these ideas in view of 
recent findings. In the case of La 2 - :I ;Sr ;r Cu04 the vHs is 
known to be crossed in the deeply overdoped region at 
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14j . A similar situation occurs 
with Bi 2 Sr 2 CuOe[15j and deeply overdoped Tl 2 Ba 2 CuOg 
also lies close to the saddle-point vHs[l6|. We will re- 
turn to these systems later. They are all single Cu0 2 
layer compounds. Where there are two Cu0 2 layers per 
unit cell the weak electronic coupling between the layers 
lifts the degeneracy of the electronic states in the lay- 
ers causing split antibonding and bonding bands. The 
splitting is maximal near (n, 0) and is about 100 meV 
there. Surprisingly it has recently been shown by ARPES 
measurements that for Bi 2 Sr 2 CaCu 2 0s the antibonding 
vHs is crossed in the deeply overdoped region around 
p ~ 0.2250. One has therefore to consider the hypothe- 
sis that the location of a vHs around p ~ 0.23 is general 
and possibly plays a central role in defining the generic 
phase curve T c (p) (along with the pseudogap at p cr u)- 
This is given all the more weight by the fact that both 
single- and double-layer cuprates cross the singularity at 
about the same point in the phase diagram in spite of 
the split band in the latter case. 

We have explored this hypothesis using a rigid ARPES- 
derived E(k) dispersion for Bi 2 Sr 2 CaCu 2 0s expressed 
in terms of a six-parameter tight-binding fit to the 
dispersion Q. The doping dependence of the resultant 
density of states at the Fermi level, N(Ep), is shown in 
Fig. QJa). Elsewhere we showed that the entropy and 
superfluid density calculated from this dispersion agrees 
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with the observed magnitude, temperature and doping 
dependence of these parameters 17f . Here we ignore the 
pseudogap in underdoped samples and this is a feature 
that would need to be added in a fuller treatment (though 
the general conclusions would remain unchanged). 

In the superconducting state we employ a d-wave gap 
of the form A^. = ^A g-^ where g-^ = cos k x — cos k y . The 
dispersion in the presence of the superconducting gap is 

given by E^ - 



e^. + A^, and Ao (T) is determined from 
the self-consistent weak-coupling BCS gap equation [l8j 
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For a justification of this approach see Ref. [171 ] . We 
adopt a pairing potential of the form = Vg-^g-^. 

The summation in Eqn. [5] extends over all states up to 
a constant energy cut-off w c =125meV. We consider two 
cases: (i) the pairing amplitude, V, is constant, chosen 
such that T c , max takes the observed value; and (ii) for 
each p- value V is selected such that T c (p) follows the 
experimentally-observed, approximately parabolic, phase 
curve. Fig.QJa) shows T c (p) plotted as a function of hole 
concentration for the two cases. 

Turning first to the T c (p) curve for a constant V (up- 
triangles), it is evident that if the pairing interaction or 
cut-off energy is fixed then the phase curve (a) is more 
narrow than that which is observed, (b) maximises at 
the location of the vHs in the heavily overdoped region 
(not at optimal doping) and (c) exhibits a second peak at 
the bonding-band vHs. All three difficulties are averted, 
and the peak broadened and shifted back to the observed 
optimal doping, only if either the pairing interaction or 
the cut-off decreases rapidly with doping. This is a robust 
result independent of the particular details that follow. 
We illustrate this by the second case explored, as follows. 

The second T c (p) phase curve (squares) follows the em- 
pirical, approximately parabolic, phase curve Q 



82.6(p-0.16) 2 



(3) 



Using the gap equation, we have calculated values of V 
that reproduce these T c values. The resultant values of 
V are shown by the circles plotted as a function of p 
in Fig. Qlb). They descend rapidly towards zero with 
increasing doping. If, alternatively, V is held constant 
and the energy cut-off, uj c , is varied, essentially the same 
result is obtained - a rapidly descending value that van- 
ishes near p rj 0.3. Such a strongly dependent interaction 
would not usually be associated with the electron-phonon 
interaction where the phonon energy scale is only a weak 
function of doping. Moreover, for a phonon mechanism 
involving motion of atoms lying outside of the CuC>2 
plane, this function is likely to be non-universal due to 
the different bond lengths and substitutionary doping 
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FIG. 1: (Color online) (a) the doping dependence of the DOS 
N(Ef) and of T c , as observed (squares) and as calculated 
(triangles) assuming a constant pairing interaction, (b) the 
doping dependence of the pairing amplitude, V, and of the 
exchange energy J, the pseudogap energy E g , the NMR wipe- 
out line kBT wo and Neel line fcsTjv. 



mechanisms in each high-T c superconductor. These re- 
sults are, rather, indicative of a magnetic, or magnetically 
enhanced, mechanism. 

Further, the pairing pairing potential, V, is clearly 
large and grows with underdoping towards the magni- 
tude of J the exchange interaction. For comparison, we 
plot in Fig. [TJb) the magnitude of J determined from 
two-magnon Raman scattering (l9j. (where J is taken, 
as usual, as 1/3 the frequency of two-magnon scatter- 
ing peak). The magnitude and doping dependence of V 
is very similar to that of J, suggesting a close relation- 
ship between these. Also plotted (open circles) is V when 
c<j c =100meV. A similar rapid rail is found showing that 
the choice of lu c is not too critical. 

This rapid fall in energy scales with doping is also 
reflected in serveral other energy scales also shown in 
Fig- QIb). These are the pseudogap energy scale, E g and 
line T wo /kB where T wo is the temperature where NMR 
intensity wipeout effects are observed, indicating the on- 
set of inhomogeneous spin and charge distribution (2p|. 
2l|. T wo is for (Y,Ca)Ba2Cu 3 O7_*0. It is also in 
this region that the 4x4 checkerboard structure is ob- 
served in scanning tunneling spectroscopy [2^. These 
lines all expand out from the antiferromagnetic phase 
curve, Tm (v) [23| . like ripples of remanent magnetic ef- 
fects suggesting a common magnetic origin for these. 
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Here the doping is estimated from the parabolic phase 
curve which we know to be approximate only. In fact 
V is very linear in Ep — E v h s1 as shown in the inset 
to Fig. [2(a). This suggests that the overall phase curve 
T c (p) is indeed governed by the proximate vHs combined 
with a rapidly declining bosonic energy scale. The value 
of uj c or V need not vanish at p ss 0.27. Eventually the 
superconducting energy gap will fall below the pairbreak- 
ing scattering rate and T c will be reduced to zero [11] even 
if uj c or V are not quite zero. 

If the overall phase diagram in the above-noted single- 
layer cuprates and in Bi-2212, is controlled by the prox- 
imate vHs then one might expect a similar vHs cross- 
ing in strongly overdoped Yo.8Cao.2Ba 2 Cu307_,5. Here 
we now present evidence for this from the T-dependent 
89 Y Knight shift measurements reported previously by 
Williams et aL[25( and shown in Fig. [2] The Knight shift 
is linearly related to the spin susceptibility, Xs> as follows 



200 



Q_ 
Q_ 



a: 



100 



-100 - 



-200 




100 200 300 

Temperature (K) 



400 



FIG. 2: The 89 Y Knight shift for Yo.8Cao.2Ba2Cu 3 7 -,5 at the 
indicated p-values. Solid curves are the fits using a rigid two- 
band dispersion, as described in the text. The inset shows the 
deduced values of the pseudogap energy, E g , and of Ef—E v h s 
showing the vHs is crossed in the deeply overdoped region. 
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K(T) = a X s(T) + a c 



(4) 



where a c is the T- and p- independent chemical shift. For 
a Fermi liquid the spin susceptibility is 



(5) 



Xs = 2mb / (df/dE)N(E)dE. 



Here /ib is the Bohr magneton and f(E) is the Fermi 
function. We adopt our previous model which we used 
to model the entropy and superfluid density [17|. Using 
a tight-binding two-band fit (provided by the authors of 
Ref. [3]) to the electronic dispersion for YBa2Cu30y_5 
and assuming a rigid band structure with doping we have 
calculated Xs an d fitted to the Knight shift data shown 
in Fig. [2] As before 17| , we adopt a Fermi-arc model 
of the pseudogap and ignore all states in the integral 
with energy below the fc-dependent gap value E g (8) . The 
fits are shown by the solid curves in the figure and the 
deduced values of E g (—E g (9 = 0) and Ep — E V Hs are 
plotted in the inset as a function of doping. 

Just like S el /T, where S el is the electronic entropy, 
the Knight shift is observed to exhibit a fanning out be- 
havior as T is decreased, with underdoped samples de- 
creasing (due to the pseudogap) and overdoped samples 
increasing (due to the proximate vHs)[3, In fact 

it has previously been shown that Xs ~ awS el /T over 
a broad range of T and p for a number of cuprates, 
with aw being equal to the Wilson ratio for nearly- 
free electrons [3|. With increasing doping both Xs(p) and 
S el (p)/T should rise steadily to a peak at the vHs then 
fall, as has been observed in La 2 _ 2 ,Sr 2 ,CuO40. The ap- 
parently classical increase in S el (p) of 1 fcs per added 
holei would appear simply to arise from the increased 
DOS on approaching the vHs. The inset to Fig. [2] shows 
that for Y .8Cao.2Ba2Cu307_5 the Fermi level will reach 
the vHs in the deeply overdoped region, just like the other 
cuprates. 



We therefore believe this situation to be universal in 
the cuprates. It raises the issue as to whether the phase 
diagram is truly a universal function of hole concentra- 
tion, p, or whether it is possibly merely a universal func- 
tion of the distance away from the vHs i.e. Ep — E v h s - 
From our point of view this would be a significant con- 
cession as we have long argued for a universal T-p phase 
diagram 0, . But, it is not straightforward to deter- 
mine doping levels except in the case of La2- a; Sr 3 ;Cu04 
and La2- K Sr :E CaCu206 which both follow the universal 
phase curve [4(. To these we can add the fully deoxy- 
genated material Yi_ a; Ca a ;Ba2Cu306 where p = x/2, ex- 
actly. Here the universality of the phase diagram seems 
to be preserved [1| and the careful use of bond valence 
sums also seems to confirm this for Yi_ ;r Ca a ;Ba2Cu307_ ( 5 
when 5 < l[3|. 

Turning to other cuprates, very precise measure- 
ments of the cation content of oxygen-stoichiometric 
Tlo. 5 Pbo. 5 Sr2Ca 1 _ a; Y 2; Cu20 7 , Tl .5Pbo. 5 Sr 2 Ca 2 Cu 3 9 
and Ti2Ba2Ca2Cu30io also confirm the basic phase di- 
agram with optimal doping at p « O.160. For other 
cuprates we have usually estimated the doping state 
either from the room-temperature thermopower [2J or 
from the parabolic phase behavior given by Eqn. [3] 
HgBa2Cu04+5 is also found to be consistent with both 
of these correlations [27j . But we have recently found 
that the overall doping dependence of the ther mop ower, 
Q(T,p), is governed by the approach to the vHs[28], with 
the change in sign of Q(T,p) at low temperature occur- 
ring precisely at the vHs, p = p V Hs- Thus Q(T,p) is 
primarily governed by (Ep — E v h s )- The thermopower 
is therefore only a universal function of p if (Ep — E v h s ) 
is a universal function of p. This can be tested by de- 
termining the absolute doping state from the area of the 
Fermi surface. We will call this pfs- 

Consider a few examples. (i) Firstly, it has al- 
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ready been shown that La2- x Sr x Cu04 appears to fol- 
low the Luttinger theorem and pps ~ x across the 
phase diagram[lj|; (ii) ARPES measurements on fully- 
oxygenated YBa2Cu3 0g.993, which we have previously 
identified as having a doping state of p=0.19[24j, has a 
Fermi surface area of 60% of the Brillouin zone corre- 
sponding to p=0.20; and (iii) the Fermi surface of Tl- 
2201 has been measured by angular magneto resistance 
oscillations [2§| and for a sample with T C =30K pps wa s 
found to be 0.24 ±0.02 holes/Cu. Using Eqn. E| its T c 
value implies p=0.25. These successes looks very promis- 
ing. However, in the case of Bi-2212 the doping state 
is a little higher than estimated from thermopower or 
Eqn. E]and for Bi-2201 it is seriously higher^. We thus 
conclude that it is likely that the phase diagram is prob- 
ably a universal function of (Ep — E v h s ) rather than of 
p (though it is often also a universal function of p) . 

Finally, we note that close to a vHs there should be 
a tendency to structural change or local symmetry re- 
duction so as, for example, to split the vHs and thus to 
lower the electronic energy. In this regard it is interest- 
ing that most, perhaps all, of the above mentioned HTS 
systems tend to be unstable in the heavily overdoped re- 
gion. Firstly, the vHs is evidently close to the limits of 
overdoping of Bi-2212 and Y,Ca-123. Our experience is 
that these systems tend to decompose near these limits. 
Thus we typically have to oxygenate at high oxygen pres- 
sures using rather low temperatures (< 350°C) to avoid 
decomposition. For a long time La-214 has been known 
to phase separate in this overdoped region[3l| and only 
by quenching from high temperature synthesis and reoxy- 
genating (to stoichiometric 0=4) at substantially lower 
temperatures can this phase separation be avoided [32j|. 
These effects are most noticeable in polycrystalline sam- 
ples where there is a large surface to volume ratio. Single 
crystals enjoy a metastable state much longer or to higher 
temperatures. Oxygenation of Tl-2201 into the heavily 
overdoped region results in precipitation of TI2O3 on the 
surface of single crystals or in grain boundaries of poly- 
crystalline samples. Our attempts to overdope Tl-2201 or 
Tlo.5Pbo.5Sr2CaCu207 into this region by Cd substitu- 
tion for Tl (generally successful in other circumstances) 
failed. And then Bi-2201 undergoes a change in elec- 
tronic state in this same region. These various instabil- 
ities have always been a puzzle as has the fact that so 
many cuprates, as prepared, reside close to optimal dop- 
ing. The proximate vHs could be the common cause. 

In summary, we have used a rigid dispersion for 
Bi2Sr2CaCu20g+5 and computed the phase curve T c (p) 
showing that the only way in which the observed phase 
curve may be recovered is if the pairing interaction is 
a rapidly falling function of doping. We find that this 
closely follows the doping dependence of the exchange in- 
teraction, J. We also argue that all HTS cuprates exhibit 
a van Hove singularity in the deeply overdoped region 
and this probably underlies the universal phase diagram 



observed in these systems. 
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